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ABSTRACT 

The census of heavy elements (metals) produced by all stars through cosmic times up to present-day is limited to < 50%; of these 
only half are still found within their parent galaxy. The majority of metals is expelled from galaxies into the circumgalactic (or 
even more distant, intergalactic) space by powerful galactic winds, leaving unpleasant uncertainty on the amount, thermal properties 
and distribution of these key chemical species. These dispersed metals unavoidably absorb soft X-ray photons from distant sources. 
We show that their integrated contribution can be detected in the form of increasing X-ray absorption with distance, for all kinds 
of high-energy cosmic sources. Based on extensive cosmological simulations, we assess that ~ 10% of all cosmic metals reside in 
the intergalactic medium. Most of the X-ray absorption arises instead from a few discrete structures along the line of sight. These 
extended structures, possibly pin-pointing galaxy groups, contain million degree, metal-enriched gas, 100-1,000 times denser than 
the cosmic mean. An additional ~ 10% of cosmic metals could reside in this phase. 

Key words, intergalactic medium. X-rays: general, gamma-ray burst: general, galaxies: active, galaxies: halos 


1. Introduction 

Elements heavier than boron (in short, “metals”) are produced 
only by stars and supemovae. Given the star formation history 
of the universe (Madau & Dickinson 2014) one can compute the 
amount of metals produced and then dispersed in the surround¬ 
ings by supernova (SN) explosions (core-collapse SNe provide 
~ 80% of the metals by mass) and stellar winds from asymptotic 
giant branch stars. The relative distribution of metals among dif¬ 
ferent galactic components (galaxy, halo, intergalactic medium 
[IGM]) is a sensitive diagnostic of such complex phenomena 
as star formation, galactic outflows, and galaxy interactions that 
might strip metals from their formation sites. 

Metal reservoirs within galaxies are stars, interstellar 
medium, and dust (Fukugita, Hogan & Peebles 1998). At the 
present epoch (redshift z = 0) the fraction of metals produced 
by star-forming galaxies and retained within them is ~ 20 - 25% 
(Peeples et al. 2014; Shull, Smith & Danforth 2012). This calls 
for an efficient transport mechanism that pours metals from 
galaxies into the circumgalactic medium (a few hundreds kpc 
size environment also often referred to as the galactic halo) or 
the “true” IGM. These galactic outflows can be powered ei¬ 
ther by SN-driven galactic winds (Oppenheimer & Dave 2006; 
Kobayashi, Springel & White 2007) or by energy injection from 
an active galactic nucleus (Planelles et al. 2014; Barai et al. 
2013). This energy deposition may also strongly suppress the 
galaxy star formation activity. 

Galactic haloes are characterised by a wide range of densi¬ 
ties and temperatures. Metals in this medium are mainly detected 
through optical/UV absorption lines of intervening absorbers 
along the line of sight (LOS) to bright sources, typically quasars. 
A considerable observational effort has been put into these stud¬ 


ies. The most prominent absorption features correspond to high- 
ionization species (indicating the presence of warm gas with 
T ~ 10^ K) like OVI, CIV, SiIV (Tumlinson et al. 2011; 
Danforth et al. 2014; Werk et al. 2014), and dust (Menard et 
al. 2010; Menard & Fukugita 2012). These metals, which can be 
traced up to galactocentric distances of ~ 150 - 500 kpc, add 
~ 20 - 30% to the total metal census, still leaving about 50% of 
the total metal budget undetected. This is known as the miss¬ 
ing metals problem (Pettini et al. 1999; Ferrara, Scannapieco 
& Bergeron 2005; Bouche, Lehnert, Peroux 2005; Pagel 2008; 
Ford et al. 2014). 

The missing metals problem looks somewhat different at dif¬ 
ferent redshifts. The metal census at z ~ 2 has been analysed in a 
series of papers by Bouche et al. (2005, 2006, 2007). Currently 
known galaxy populations at z ~ 2 can account for 30 - 60% 
of the expected metals. Absorption line studies trace the metal 
content in lower density gas phases. Metals in the Lya forest 
(logAni < 17) account for ~ 15 - 30% of the metals. Lyman 
limit systems (log Ari > 17, including Damped Lya systems), 
provide ~ 5 - 15% (Lehner et al. 2014; Bouche et al. 2007) 
and at least 10% is in the form of highly ionised metals in the 
circumgalactic gas of galaxies (Lehner et al. 2014). Summing 
up these contributions we conclude that 50% of metals lies out¬ 
side galaxies. Based on these estimates, one can conclude that 
the metal census at this redshift might be almost complete. On 
the contrary, in the local universe, Peeples et al. (2014) showed 
that - 50-55% of metals can be accounted for within galaxies, 
with an error budget of - +25%. The remaining metals should 
be dispersed outside galaxies on a larger scale. 

Thus we are forced to make the hypothesis that the missing 
metals should be found at even larger distances, in the IGM or 
in the gas in between galaxies in groups and clusters. Among 
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Sample 

Instrument 

Detections (UL) 

Redshift range 

Refs. 

GRB-C1012 

Swift 

107(11) 

0.03-8.2 

1,2 

AGN-E13 

XMM-Newton 

24 (34) 

0.45-4.72 

3 

AGN-Cll 

XMM-Newton 

40 (245) 

0.02-2.39 

4 

AGN-Y09 

XMM-Newton 

32(141) 

0.22-5.41 

5 

AGN-hz 

XMM Sl Chandra 

1(15) 

3.63-6.28 

6,7 

AGN-P05 

XMM-Newton 

19 (44) 

0.01-6.28 

8 


Table 1. Summary of X-ray source samples used in the text. 


References: 1: Campana et al. (2010); 2: Campana et al. (2012); 3: Eitan & Behar (2013); 4: Corral et al. (2011); 5: Young, Elvis & Risaliti (2009); 
6: Vignali et al. (2003); 7: Saez et al. (2011); 8: Page et al. (2005). Superposition of data were checked and older data were removed. 


these possibilities, we can safely exclude only clusters, as their 
metal content is relatively easy to measure and largely insuf¬ 
ficient to close the gap (Simionescu et al. 2011; Werner et al. 
2013). Observations of galaxy groups, although much more dif¬ 
ficult, tentatively indicate a larger metal content than the sum of 
single galaxies (Sun et al. 2009; Humphrey et al. 2012). Finally, 
little can be said about the amount of metals in the IGM. 

Metals absorb soft X-ray photons. This “metal fog” should 
manifest itself as an increase in X-ray absorption with distance 
for all kinds of cosmic high-energy sources. A first hint of the 
existence of this metal fog came from gamma-ray bursts (GRB) 
observations (Campana et al. 2010; Behar et al. 2011). These 
sources, spanning a very large redshift range, are the primary 
probe for this effect. Despite the intrinsic absorption by gas 
within the GRB host galaxy itself (see Fig. 1), a minimum ab¬ 
sorbing column density envelope enclosing all GRBs and in¬ 
creasing with redshift exists. If it arises from a cosmological 
gas distribution, this minimum cumulative absorption could af¬ 
fect all kinds of high-energy sources. Indeed, this effect persists 
when GRB data is combined with active galactic nuclei (AGN) 
samples(Starling et al. 2013; Eitan & Behar 2013). At variance 
with other techniques (e.g. UV absorption line studies), X-ray 
absorption is less sensitive to the ionisation state and thermo¬ 
dynamic phase (gaseous vs. solid, i.e. dust) of the absorbing 
medium, allowing for, in principle, an almost complete metal 
census. 

In this paper we explore the minimum X-ray absorbing col¬ 
umn density based on samples of GRBs and AGNs discussed 
in Section 2 (particular care is taken to discuss a few out¬ 
liers present in these samples and discussed in Appendix A). In 
Section 3 we work out the lower envelope of the column den¬ 
sity distribution as well as the mean absorbing value. We then 
investigate how this absorption arises by means of cosmological 
simulations. In Section 4 we briefly describe our simulations. In 
Section 5 we describe in detail how we evaluate the absorbing 
cross section of the intervening material. The evaluation of the 
fraction of metals in the IGM and in the cosmic structures is car¬ 
ried out in Section 6. Discussion and conclusions are drawn in 
Section 7. 

2, Sample of X-ray sources 

Our sample of absorbed X-ray sources consists of GRB and 
AGN. We took data from the literature without further interven¬ 
tion (expect for a few outliers, see the Appendix). We selected 
data in order to have a sufficiently large number of counts to 
provide a reliable measure of the X-ray spectrum and of the in¬ 
trinsic absorption column density at the source redshift (A^nfe)) 
in particular. The GRB sample relies on Swift data and the AGN 
samples on XMM-Newton (see Table 1). 


2.1. Gamma-ray bursts 

We constructed the GRB sample (GRB-C1012) by merging the 
sample compiled by Campana et al. (2010) with all GRB with 
known redshifts at that time together with the additional bursts 
quoted in Campana et al. (2012), part of the complete sample 
of bright GRB (Salvaterra et al. 2012). This sample is selected 
based on a cut at the burst peak flux (BAT6) and is free of bias 
in the column density measurements, which is almost complete 
in redshift. The total number of GRB amounts to 118. We have 
an upper limit on the Nuiz) for only 11 of them, testifying that 
intrinsic absorption is a common characteristic of these sources. 
It is apparent that higher-redshift GRBs are characterised by a 
larger N^iz). Based on the Campana et al. (2010) sample, a 
Kolmogorov-Smirnov (KS) test shows the probability that the 
column density distributions of GRBs above and below z = I 
are not drawn from the same parent distribution is 0.08 per cent 
(equivalent to Gaussian 3.3 cr). This supports a cosmological in¬ 
crease of the absorbing column density (Campana et al. 2010; 
Behar et al. 2011; Starling et al. 2013). 

2.2. Active galactic nuclei 

AGN data were taken from selected bright samples of objects: 

AGN-E13: Eitan & Behar (2013) constructed a sample of 
58 high-redshift (0.45 < z < 4.72) AGN taken from the XMM- 
Newton archive and selected to have > 1, 800 counts. This en¬ 
sures an accurate determination of the intrinsic column density. 
The sample includes radio-quiet and radio-loud quasars (30 and 
28, respectively). The sample contains 24 detections and 34 up¬ 
per limits. 

AGN-Cll: Corral et al. (2011) selected AGN from the 
XMM-Newton bright survey (XBS). The XBS is composed of 
two flux limited samples selected in the medium (0.5-4.5 keV) 
and hard (4.5-7.5 keV) energy bands. The sample contains 270 
Type I (unabsorbed) AGN, 29 Type II (absorbed) AGN and five 
blazars. The distribution in redshift covers the range 0.02-2.39. 
The sample contains 40 determinations of Ah(z) and 245 upper 
limits. 

AGN-Y09: Young, Elvis & Risaliti (2009) built an AGN 
sample matching quasars from the 5^^ Sloan Digital sky survey 
release with serendipitous XMM-Newton observations. The ini¬ 
tial sample comprised 685 X-ray source detections, with 589 
radio-quiet, 62 radio-loud, and 34 broad absorption line (BAL) 
quasars. A meaningful Ah(z) can be derived for only a minor¬ 
ity of them. The final sample consists of 32 detections and 141 
upper limits. The redshift range sampled by these AGN is 0.22- 
5.41. 

AGN-hz: In order to increase the statistics at high redshift, 
we merged the small groups of sources provided by Vignali et al. 
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Fig. 1. Column densities of GRB and AGN of our sample as a function of redshift (1 + z). Black symbols are used for GRB (filled 
circles for detections and arrows for upper limits). Filled red circles and arrows identify detected (upper limits) intrinsic column 
densities within the AGN-E13 sample, which has been used to compute the mean absorbing column density. Green symbols (filled 
circles and arrows) identify (detections and upper limits, respectively) AGN of the other samples. Circled crosses indicate AGN that 
were discarded from the overall sample: magenta symbols identify z < 0.4 AGN, light blue symbols were further investigated in the 
Appendix. The continuous line marks the lower envelope of the intrinsic column densities of log Ah - 19.5 and the dashed line the 
mean of the AGN-E13 sample, as derived from the censored analysis, with log Ah - 20.3. 


(2003) and Saez et al. (2011). This is based on XMM-Newton and 
Chandra observations. In this sample, spanning 3.63 < z < 6.28, 
we have just one Ah(z) determination and 15 upper limits. 

AGN-P05: Page et al. (2005) selected a sample of 29 high- 
redshift (z > 2) AGN, including 7 radio-quiet, 16 radio-loud, and 
6 BAL quasars, and complemented it with a number of low red- 
shift AGN (39). There is a strong overlap with previous samples 
and only 59 AGN were retained (0.01 < z < 4.41). Of those, 19 
have a determination of Ah(z) and 44 only an upper limit. 

3. Lower envelope and mean of the column density 
distribution with redshift 

The total sample of objects comprises 711 sources spanning 
the redshift range 0.01-8.2. Of these, 222 yield positive detec¬ 
tions of an intrinsic column density. We approached this distri¬ 
bution with two main drivers. The first is to find out the mini¬ 
mum column density allowed by all these data. This minimum 
Ah can be derived as the envelope obtained taking a local ab¬ 
sorber and evolve it at different redshifts. In this case we are not 
interested in a fit of the different column densities, which pro¬ 


vide a very biased sample (e.g. due to different selection criteria 
and intrinsic column densities, etc.), but a lower universal en¬ 
velope applicable to all X-ray sources. It has been shown that 
a local absorbed Ah(z = 0) scales with redshift as (1 -r z)^'^ 
(Campana et al. 2014). This means the effect of an absorber 
a z = z on the X-ray spectrum of a distant source scales as 
Ah(z = 0) = Ah(z = z)/(l + The scaling index has 
been shown to slightly depend on the observing instrument with 
a small spread of +0.05 around the theoretical value of 2.40 
(Campana et al. 2014). This minimally impacts the Ah evalu¬ 
ation. To evaluate the minimum column density, we consider a 
cold absorber (i.e. with the highest cross section) with present- 
day column density Ah(z = 0) and use a value larger by a factor 
(1 +z)^'^ at any z > 0. After a careful analysis of a few outliers in 
the Ah -z plane (see above), we fit the data with the requirement 
to leave below the curve 1% of all the data (mainly provided by 
upper limits). Upper limits are treated here as detections. We find 
a value of log Ah(z = 0) ^ 19.5 (we express column densities in 
units of cm“^, see Eig. 1). 

Concerning the evaluation of the mean contribution along 
the line of sight we adopted a different approach. We stick to 
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Fig. 2. Distribution of equivalent hydrogen column densities, 
Nil, with redshift as measured for different sources: GRBs with 
measured redshift (grey dots), bright GRBs in a redshift com¬ 
plete sample (black dots), and AGN (green squares). Error 
bars are at 1 cr confidence level. The long dashed curves rep¬ 
resents the LOS with the lowest absorption, i.e. with no ab¬ 
sorbers with density contrast A > 100, and constant metallicity 
(Z ~ 0.03 Zq). This curve was obtained from cosmological sim¬ 
ulations, summing up the contribution by each cell, with an ab¬ 
sorbing column density weighted for its effective (temperature- 
dependent) value. The results for a physically motivated model 
of SN winds (momentum-driven model by Oppenheimer et al. 
2012) are shown by the solid line (median of LOS distribution) 
and corresponding shaded areas (68% and 90% of LOS, respec¬ 
tively). The short dashed line corresponds to the median LOS 
with log Ah = 20.3 + 0.1 (1 cr confidence level) obtained consid¬ 
ering constant metallicity (Z = 0.3 + O.IZ©) in absorbers with 
100 < A < 1000. 


a well-selected and characterised sample and with high statis¬ 
tics (> 1,800 counts per spectrum), which is the AGN-E13 
sample (Eitan & Behar 2013). In this sample only 43% of the 
quasars show significant absorption. The authors tried to find 
which physical parameters of the quasars might drive the ob¬ 
served absorption, e.g. time (i.e. evolution), radio luminosity, or 
the X-ray luminosity. The most likely interpretation is in terms 
of an increasing absorption along the line of sight, possible with 
minimal contamination from intrinsic absorption in the AGN 
themselves. We note in fact that there are no highly extinguished 
objects in the sample, hinting at the lack of Type II, intrinsically 
heavily absorbed, AGN. 

We made minimal changes to the AGN-E13 sample, simply 
getting rid of the two outliers above. The sample then consists of 
56 objects. More than half of the Ah estimates are upper limits 
(31, 55%). This large fraction of upper limits calls for censored 
statistics. It is not easy to perform a linear regression test, keep¬ 
ing fixed the value of the slope (to 2.4). Lor this reason we recti¬ 
fied all the Ah estimates by dividing them for the term (1 -r z)^-^. 
We then approached the new distribution, searching for the best 
(censored) constant fit. We used programs of the ASURV suite 


(Isobe et al. 1986). In ASURV we made use of the Kaplan-Meier 
estimator. This censored fit results in a mean column density ab¬ 
sorption of log Ah(z = 0) = 20.3 + 0.1 (1 cr confidence level). 

If we restrict ourselves to a smaller sample, obtained by re¬ 
quiring at least 5,000 counts, we retain 24 objects, with only 5 
upper limits (21%). In this case where the impact of upper lim¬ 
its is much lower, we obtained a minimum column density of 
^19.3 and a mean column density absorption of ^20.4. 


4. Cosmological simulations 

We have used the simulation presented in Pallottini et al. (2013), 
whose main features are recalled below. The simulation is per¬ 
formed using the publicly available adaptive mesh refinement 
code RAMSES (Teyssier 2002), which can be described as a fully 
threaded tree data structure where the hydrodynamical adaptive 
mesh refinement scheme is coupled with a particle mesh N-body 
solver employing a cloud-in-cell interpolation scheme to solve 
the Poisson equation. 

Note that heating-cooling processes in RAMSES are already 
implemented (Theuns et al. 1998; Aubert et al. 2008), including 
the treatment of an external UV ionising background (Haardt & 
Madau 1996); we have not included star formation, metal enrich¬ 
ment, nor internal radiative feedback for the following reasons. 
The treatment of winds, metal transport, and mixing remains an 
extremely difficult challenge for all current simulations. Hence 
we have decided to use a simulation containing only physics that 
we can reasonably keep under control (gravitation, structure for¬ 
mation shocks and radiative processes, adiabatic cooling), and 
combined it with different prescriptions for the metal spreading 
taken from specific studies or proposals by other leading groups. 
This appears reasonable as our aim is to bracket the uncertain¬ 
ties of the problem rather than constraining different feedback 
and enrichment models. 

The simulation box size of 100 h~^ (h is the normalised 
Hubble constant in units of 100 km s“^ Mpc“^) comoving Mpc 
is evolved from z = 100 to z = 0 using 256^ dark matter (DM) 
particles and a corresponding number of base baryonic grid; we 
allow six additional levels of refinement adaptively chosen with 
a Lagrangian mass threshold-based criterion. This yields a mass 
resolution of 1.65 x 10^^ Qdm h~^ Mq for the DM and a formal 
spatial resolution of 6h~^ kpc for baryons. We have assumed 
a WMAP7 cosmology (Larson et al. 2011) with Qa = 0.727, 
Qdm = 0.228, Qb = 0.045, n = 0.967, erg = 0.811, and 
h = 0.704. 

Simulation results, in particular the baryon thermodynamic, 
are consistent with numerical studies (Rasera et al. 2006; Peeples 
et al. 2010). Lor example, Peeples et al. (2010) simulated a 
12.5h~^ Mpc comoving box with a high (2 x 288^ particles) 
and low (2 x 144^ particles) resolution using the smooth par¬ 
ticle hydrodynamic (SPH) code GADGET-2 (Springel 2005). At 
z = 3, Peeples et al. found an equation of state very similar to that 
obtained here. Additionally, when restricting to the quasi-linear 
regime, baryons evolution is in accordance with semi-analytical 
models (Miralda-Escude, Haehnelt & Rees 2000; Choudhury, 
Padmanabhan & Srianand 2001). 

Lor the present work we have extracted various LOS from the 
light cone going from z = 0 to z = 2. As in Pallottini et al. (2013), 
the LOS are randomly selected from contiguous snapshots of the 
simulation. Lor this work a suitable convergence of the results is 
obtained with A = 100 LOS. 
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Fig. 3. Individual contribution of cells in the At-Tt plane for the reference momentum-driven wind model along three different LOS: 
Left panel a) refers to the LOS with the lowest absorption. The optical depth along this LOS in our simulation is T(0.5keV) = 0.06. 
Central panel b) indicates the median LOS corresponding to the solid, black line in Fig. 2, the median optical depth is T(0.5keV) = 
0.21. The right panel c) refers to the LOS defining the top envelope of simulated curves in Fig. 2. In this case the optical depth is 
T(0.5keV) = 0.88. The colour scale represents the value of t/ computed at 0.5 keV. The LOS with the lower absorption lacks of 
cells with Ai > 100. Instead these cells dominate the X-ray optical depth in the median and top LOS, where most (~ 60%) of the 
absorption is provided by two systems with A/ > 300 and T/ > 10^ K (dark red points). 


5. Cosmological absorption 

We compute the X-ray optical depth for photons seen by the 
observers at the energy E by summing up the contribution of 
each cell along the LOS up to the source redshift. The i-th con¬ 
tribution at the energy E' = £’(1 -r z/) is given by rfE') = 
^iNuiCriE', where A^h,, Tf are the gas column density and 

temperature of the i-th cell and are provided directly by the sim¬ 
ulation. The ionisation parameter is ^ ^TrEcxRBMei^ where 
the integrated cosmic X-ray background spectrum is Ecxrb ~ 
2.9 X 10“^ erg cm“^ s"^ sr“^ in the energy band 0.005-300 keV 
(Deluca & Molendi 2004) and the election density of the i-th 
cell, ne-, is a factor -1.2 higher than the hydrogen density in 
the cell. Finally, we computed the photo-electric cross-section 
cr using the XSPEC model ABSORI (Done et al. 1992) for the 
temperature £/ and ionisation parameter , assuming the solar 
composition on a parameter grid. We assigned the metallicity of 
the i-th cell, Z/, according to different assumption and physical 
motivated models (see next Section). Therefore, we consider the 
temperature and ionisation state for each cell, depending on the 
density, and we work out the appropriate cross section. This is 
computed at 20 different (z = 0) energies logarithmically spaced 
in the 0.3-10 keV energy range. All these contributions are then 
summed up at the given redshift to provide the total absorbing 
column density along the line of sight. To compute the intrinsic 
Nil to be compared with observational data, we fit the resulting 
cross section with a cold absorber with ^ = 0 at the redshift of 
the background source. 


6. Metal enrichment models 

The simulation described in Sect. 4 does not directly provide the 
metallicity of the gas that is needed to compute the X-ray optical 
depth. Metals produced by supernova and during the asymptotic 
giant branch phase are expected to be spread in the circumgalac- 
tic and intergalactic medium by galactic outflows (Oppenheimer 
& Dave 2008; Kobayashi, Springel & White 2007). In the fol¬ 


lowing we will describe how we assign metallicity to each simu¬ 
lated cell based on its over-density (A - Z prescription) and how 
our results depend on the different choices. In particular, we pro¬ 
vide three different descriptions for the metal enrichment of in¬ 
creasing complexity, starting from back-of-the-envelope calcu¬ 
lations to physically motivated models. The estimates we derive 
are consistent. 

6.1. Simple estimates 

The photo-electric optical depth (opacity) of a homogeneous 
IGM of constant metallicity Z (in solar units) at observed pho¬ 
ton energy E and temperature T can be expressed as (Behar et 
al. 2011; Starling et al. 2013): 

dt' 

tigm{E,z.Z)= nii(z')cr(E,z\Z,T)c — dz'. (1) 
Jo dz' 

The absorption cross section (t(£, z\ Z, T) for a solar mixture of 
elements scales linearly with Z. For a homogeneous and isother¬ 
mal IGM one can write the above equation as 

t,gm{E, z, Z) = ^ { z(t{E, z',Zo, T) 

LIq Jo 


V(1 + 

where no = 1.7 x 10“^ cm“^ is the present-day mean IGM den¬ 
sity, Ho = 71 km s“^ Mpc“^ is the Hubble constant, Qm = 0.27 
and Qa = 0.73 are, respectively, the present-day matter and dark 
energy fractions of the critical energy density of the Universe, 
and c is the speed of light. We selected an energy ofE = 0.5 keV 
just below the oxygen edge. On the right hand side of Eq. 2, r 
can be expressed as r = To obtain the IGM metallicity, we 
can consider our low column density envelope logA^n = 19.5. 
In this case, however, the appropriate cross section to use on the 
right-hand side of Eq. 2 is the one pertaining to a cold absorb¬ 
ing material since all the column densities were obtained in this 
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way, i.e. cr = cr(E = 0.5 keV, z', Z©, T = 0). We can then select 
a single temperature for the IGM and work out the cross section 
dependence on redshift using the ABSORI model within XSPEC, 
as described above (this dependence is not trivial and cannot be 
expressed as a simple power law). Based on our minimum ab¬ 
sorbing column density, working out the integral (from z = 0 
to up to z = 2, it being negligible above this value) and further 
assuming an IGM T = 10^ K or T = 10^ K, we derive a mean 
metallicity of Z = 0.043 Z© or 0.025 Z©, respectively. For a neu¬ 
tral (cold) IGM the cross section on the right- and left-hand sides 
of the equation can be simplified, and we derive 0.014 Z©. 

6.2. Simple models 

To bracket all possible enriched scenarios, we consider two lim¬ 
iting cases: (i) constant metallicity, i.e. Z^/Z© - Za - const, in 
all cells, and (ii) a no-wind model. In the latter case, log(Z/) = 
\og{Z}y) -h log(A/), with Z}y =const. For both scenarios we com¬ 
pute the minimum envelope curve by considering the LOS with 
the lowest absorption (we checked that we obtain consistent re¬ 
sults using the lowest three LOS). The values of our free pa¬ 
rameters (Za and Ziy, respectively) are then obtained by requiring 
that only 1% of the data points fall below the simulated curve. 
We find: (i) Za = 0.030Z© ; and (ii) Z}j = 0.0018 Z©, the two 
curves being almost identical in Fig. 2. From this we estimate 
the metal cosmic density by eq. (1) of Ferrara et al. (2005) to be 
qigm _ 1 77 X 10“^ and 1.17 x 10“^, respectively. In spite of 
the very different (extreme) assumptions, the two models result 
in similar metal column and cosmic densities. 

In a similar manner, we compute the median of the simulated 
LOS distribution for the constant metallicity case. The signal is 
in this case dominated by over-densities in the range 100-1000. 
To match the observed median found in Sect. 3, a mean metal¬ 
licity of Z = 0.3 + 0.1 Z© associated with these over-densities is 
required. 

It is meaningful to compare the above numbers with the 
comic density of metals produced by SNe by z = 1 where the 
peak of the cosmic star formation rate is expected. Therefore, 
we integrate the star formation history (SFH) as computed by 
Li (2008) from high redshift (z ~ 15, well enough to assure 
convergence) down to z = 1, assuming that the specific frac¬ 
tion of heavy elements produced by stars is 0.03 (Peeples et 
al. 2014). Uncertainties on the integrated SFH should be minor 
since it rapidly decreases at z > 2 and Pop III gives a negligi¬ 
ble contribution for z < 4 (Tomatore et al. 2007; Pallottini et 
al. 2014). Uncertainties in the nucleosynthetic yields are instead 
larger (and not easy to quantify, e.g. Wilkins et al. 2008) and we 
decided to stick to the Peeples et al. (2014) estimate (also for 
comparison purposes). We obtained Qz(z = 1) ~ 1.4 x 10“^. 
Thus, we find that 12% and 8% of the metals produced by SNe 
by z = 1 are expelled far from galaxies in regions with A < 100 
for model (i) and (ii), respectively. 

6.3. Physically-motivated enrichment models 

We assign the metallicity to each cell in our simulation by adopt¬ 
ing the distributions of metallicities as a function of both over¬ 
density and redshift predicted by three different wind models 
of Oppenheimer et al. (2012): (a) a constant wind (CW) model 
where the ejection velocity is v^ind = 680 km s“^ and mass load¬ 
ing factor of T] = 2 for all galaxies; (b) a slow wind (SW) model 
in which the ejection velocity is half of the previous velocity; (c) 
a momentum-conserved (VZW) model with Vwind = a//l - 1 


in which fi is the luminosity in units of the Eddington luminos¬ 
ity required to expel gas from the galaxy potential, t] = crolo-g, 
and (To = 150 km s“^ , and o-g is the galaxy’s internal veloc¬ 
ity dispersion (see Oppenheimer et al. 2012 and Oppenheimer 
& Dave 2008 for model details). Oppenheimer et al. (2012) 
show that in spite of the significant difiTerences among the metal 
enrichment prescriptions, all wind models produce metal-line 
statistics within a factor of two of existing HST/COS observa¬ 
tions (Danforth et al. 2014; Liang & Chen 2014; Stocke et al. 
2013; Tilton et al. 2012; Prochaska et al. 2011), with the physi¬ 
cally motivated VZW model preferred by the data. In the follow¬ 
ing we will refer to the VZW model as our reference prescrip¬ 
tion. 

For each wind model we compute the expected Vh - z distri¬ 
bution following the procedure described above. Note that since 
the metallicity is fixed, the model has no free parameters. In Fig. 
2 we show the results for our reference prescription. In the fol¬ 
lowing we will discuss the main results obtained with this metal 
enrichment recipe and the differences with the alternative wind 
models. 

In the reference prescription (momentum-driven winds) 
most of the absorption comes from regions with 10 < A < 100 
enriched to Z ^ 0.030 Z© (see Fig. 3). They can be naturally 
identified with the genuine IGM. These regions contain 6% of 
the metals produced by SNe by z = 1; an additional 4% has been 
transported into regions of even lower over-density. 

The minimum envelope of the Vh - z distribution of our ref¬ 
erence prescription reproduces very well the observations with 
< 1% of data falling below the lowest simulated curve. The 
metal cosmic density derived from this LOS is Qz(z = 1) = 
1.16 X 10“^, i.e. ~ 10% of the metals produced by SNe by z = 1. 
This is not true for alternative models: ~ 3% of data are found 
to be outliers in the CW model, while this fraction reduces to 
~ 1% in the SW model. This reflects the fact that metal enrich¬ 
ment of the regions at the lowest over-densities is much faster 
and efficient in the CW model with respect to the VZW model. 
In particular, the CW model predicts a ten times larger metallic¬ 
ity in log(A) = 0 regions than the VZW model. This is further 
confirmed by looking to the fraction of the X-ray optical depth 
provided by the cell at different over-densities shown in the right 
panel of Fig. 4 (red lines). In the CW model most of the sig¬ 
nal is provided by regions with log(A) ~ 0 while in the VZW 
model regions with 1 < log(A) < 2 dominate. The distribution 
of the SW model is intermediate between the two. This evidence 
is confirmed by looking at the distributions in terms of the gas 
temperature (right panel of Fig. 4) with cold regions dominating 
in the CW model. 

These differences become less important when looking at the 
median of the X-ray absorption distribution (black lines in Fig. 
4). Indeed, the median of the distribution is only a little sensi¬ 
tive to the assumed model. In all models the signal is dominated 
by warm-hot (6.5 < log(r/K) < 7.5), collapsed (log(A) > 2) 
regions. Most of the absorption arises from distant intervening 
sources at z > 0.5. The resulting Vh - z curves of different mod¬ 
els lie within a factor of 1.3 from our reference prescription. 

7. Discussion and conciusions 

Although sensitive to a wide range of metal ionisation states and 
species. X-ray experiments, because of limited spectral resolu¬ 
tion, lack the ability of optical/UV studies to identify individual 
absorbers along the LOS. Rather, X-ray absorption gives a cu¬ 
mulative view of metals in the z ~ 0 - 2 redshift range. We 
explore the possible existence of a minimum X-ray absorption 
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log(A) 
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Fig. 4. Fraction of the X-ray optical depth at the observed energy £" = 0.5 keV as a function of the cell over-density, temperature, 
and redshift (from top to bottom). The left panels refer to the mean of the five LOS with the minimal absorption, while the right 
panels refer to the mean of 68% of the LOS around the median. Solid, short-dashed, and long-dashed lines refer to the VZW, SW, 
and CW models, respectively. Note that the CW model slightly overpredicts the minimum envelope of the Nu - z distribution, as 
noted in the text. 


contribution due to diffuse metals. Data support the existence 
of this ’metal fog’, even if a direct measurement is prevented 
by current (and future) observational limits. In addition, if back¬ 
ground X-ray sources are projected on top of denser intervening 
regions larger column densities are retained, resulting in nearly 
an order of magnitude spread in X-ray column densities. 

A simple estimate of the metal content in the ideal case of a 
homogeneous, isothermal IGM with mean metallicity Z can be 
obtained by comparing the photo-electric optical depth with the 
observed minimum absorbing column density. For T = 10^ K 
(10^ K) we derive, by integrating in 0 < z < 2, Z = 0.043 Z© 
(Z = 0.025 Z©), corresponding to a metal cosmic density of met¬ 
als in the IGM of Of™ = 2.4 x 10“^ {Of™ = 1.4 x 10“^). For a 
neutral (cold) IGM we instead derive Z = 0.014 Z©. To provide a 
more complete description of the absorbing medium we rely on 
extensive cosmological simulations. To compute the metal col¬ 
umn density of the intervening material, we consider 100 LOS to 
distant sources through a 100 comoving Mpc Adaptive Mesh 
Refinement cosmological simulation (Pallottini et al. 2013). 

We compute the X-ray optical depth for photons seen by 
the observer at energy E as t{E) ^ + z), 

where Ah, , ct, Tf , and Zf are the gas column density, the to¬ 
tal cross section, ionisation parameter (set by the cosmic X-ray 
background), temperature, and metallicity in the i-th simulation 
cell, respectively. All these quantities are drawn from the simu¬ 
lation box, except for the metallicity Z/. As considerable uncer¬ 
tainty in numerical simulations on the transport and mixing of 
heavy elements still remains, we consider two limiting cases to 
bracket all possible enrichment scenarios: (i) constant metallic¬ 


ity, i.e. Zi/Zq = const, in all cells along the LOS; (ii) a no-wind 
model in which metals remain bound to their formation sites. 

The minimum envelope of the Ah - z distribution should 
be provided by the simulated LOS with the lowest absorption. 
Based on our simulations, this LOS does not contain any ab¬ 
sorber with over-density A/ = p//p > 100, where pi is the gas 
density in the cell and p the mean cosmic gas density. These re¬ 
gions are hot T ^ 10^“^ K and far from collapsed structures. In 
order to match the minimum envelope of the Ah - z distribution, 
we find that 8 - 12% of the metals produced by SNe by z = 1 
must have been ejected far from galaxies in the IGM. This re¬ 
sult is almost insensitive to the distribution of metals in regions 
with different over-densities. The mean metallicity of the IGM 
is Z = 0.03 Z©. More realistic prescriptions (Oppenheimer et al. 
2012) provide results within the above range (see above). Traces 
of these structures can be found through detailed X-ray absorp¬ 
tion studies (Nicastro et al. 2013). Indications of a similar metal 
content in the IGM came very recently from UV studies (Shull, 
Danforth & Tilton 2014). 

The mean of the X-ray absorption distribution based on 
AGN is even more telling. In the simulations the median LOS is 
dominated by 2-5 absorption systems with large over-densities 
(A > 300) at z ~ 0.5 - 1.2. The total column density contribution 
of these regions (at their respective redshift) is log Ah - 21.5. 
The absorbing systems have characteristic projected sizes of 
0.1-1 Mpc and are associated with dark matter halos of virial 
mass ~ (0.2 - 2) x 10^^ M©, corresponding to a virial tem¬ 
perature of ~ (3 - 15) X 10^ K. These structures are about ten 
times more massive than the Milky Way, suggesting that the ab- 
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sorption is due to the extended circumgalactic reservoir of gas 
of (small) galaxy groups rather than to single isolated galaxies. 
Given that the mean metallicity in these regions (A ~ 100-1000) 
is 0.3+0.1Z© (in agreement with current observations Humphrey 
et al. 2012 and Su 2013), we can conclude that they contribute 
an additional ~ 11 ± 4% to the metal budget. 

Direct detection of this hot, enriched material in the X-ray 
band out to the virial radius (defined such that the mean density 
of a halo of mass M is 200 times the cosmic mean) of galax¬ 
ies is challenging with current X-ray facilities (Sun et al. 2009). 
However, the presence of these reservoirs might have been iden¬ 
tified by recent optical/UV absorption studies, thus paving the 
road to deeper, dedicated studies (Mulchaey et al. 1996; Stocke 
etal. 2013). 

Acknowledgements. We acknowledges useful discussions with E. Feigelson, R 
Gastaldello, A. Moretti, & L. Sage. 
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Appendix A: Outliers 

Based on Fig. 1, we identified 11 main outliers in the distribution 
of intrinsic column densities (A^h(^))- If this is a true cosmolog¬ 
ical effect no outlier should be present. Therefore we carefully 
study these outliers in order to find an explanation for their low 
intrinsic column densities. We approached the X-ray spectra of 
these sources and tried to verify if the evaluation of the column 
density is completely satisfactory from a statistical point of view. 
In this respect we disregard AGN with an intrinsic column den¬ 
sity log Ah <19 and z < 0.3 because of a too small intervening 
length scale and, more importantly, because of the AGN influ¬ 
ence in cleaning a substantial fraction of the surrounding envi¬ 
ronment. Five AGN were discarded in this way. For the other six 
AGN, we provide a detailed description. 

A.1. QSO B0909+5312 at z=1 .377. 

QSO B0909+5312 was observed by XMM-Newton in May 2003. 
Data were retrieved from the XMM-Newton archive and repro¬ 
cessed using SAS 13.5.0. Useful XMM-Newton data were taken 
from the entire 13 ks and 17 ks exposures of pn and MOS 
detectors, respectively. Data were extracted following common 
practice and spectrally binned to 50 counts per energy chan¬ 
nel; arf and rmf files were generated using SAS tasks. Fitted 
energy ranges are 0.2-10 keV and 0.3-10 keV for the pn and 
MOS instruments, respectively. Eitan & Behar (2013) fitted the 
same data with an absorbed power-law model. They adopted a 
Galactic column density of = 1.49 x 10^^ cm“^ (Kalberla 
et al. 2005) plus and intrinsic column density at the AGN red- 
shift. Based on their spectrum they derive an upper limit of 
^^h(^) < 7 X 10^^ cm“^. We fit our data with a higher binning to 
increase statistics. We adopted the TBABS model within XSPEC 
to account for X-ray absorption (Wilms, Allen & McCray 2000) 
and with the wi Im solar abundance pattern (this is very similar 
to Asplund et al. (2009) and in the case of oxygen, the main 
driver for the absorption, is exactly the same) and the vern pho¬ 
toelectric absorption cross sections. We obtain similar results 
(A^h(^) < 5 X 10^^ cm“^, 90% confidence level for one parame¬ 
ter of interest, i.e. A^^ = 2.71) but with a reduced A^ed “ I-I^ 
(for 187 degrees of freedom, dof). Eitting the same data with a 
broken power law we obtain instead A^ed “ ^*77 (185 dof). An 
E-test indicates that the addition of the spectral break is signifi- 
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cant at the > 8 cr confidence level. The intrinsic column density 
in this case is Nuiz) < 5 x 10^^ cm“^, perfectly consistent with 
our lower envelope. 


A.2. Two AGNs in the Lockman hole at z= 1.204 and z= 1.113. 

J105239.6+572432(J1052atz = 1.113) and J105316.9+573552 
(J1053 at z = 1.204) are two Type I AGN in the Lockman 
hole with a very low absorption (Nu(z) < 6 x 10 ^^ cm“^ and 
Nu(z) < 9 X 10^^ cm“^, respectively (Corral et al. 2011). We 
took the longest XMM-Newton observation of the Lockman hole 
not affected by soft proton flare covering 100 ks. Data were ex¬ 
tracted as above and binned to 25 counts per spectral bin. 

J1052 can be fit with a simple power law, resulting in = 
1.15 (288 dof) and Nn(z) < 3 x 10 ^^ cm'^ (A^gai ^ 5 5 >< iq19 
cm“^). A better fit is obtained with a broken power law, with 
Tred “ Nu(z) < 6 X 10^^ cm“^. An F-test 

indicates that the addition of the spectral break is significant at 

> 8(T. 

J1053 falls on the dead CCD in MOSl and across dead 
columns on the pn detector. The fit with a single power law is 
good = 1.01, 129 dof) and returns A^h(^) < 2 x 10^^ cm“^. 

Both values are consistent with our lower envelope. 


A.5. 2E 1852atz=0.312. 

This Type I AGN is reported to have an intrinsic column density 
A^h(z) < 1 X 10^^ cm“^ (Corral et al. 2011). The object was ob¬ 
served by XMM-Newton on Apr 2000 in small window with the 
two MOS and in full window with the pn. Some flaring activity 
was present and filtered out reducing the exposure time to 54 ks 
and 22 ks for the MOS and pn, respectively. MOS data were fit¬ 
ted in the 0.6-10 keV range to encompass calibration uncertain¬ 
ties in small window mode. All data were binned to have at least 
25 counts per spectral bin. We fit the spectra with an absorbed 
power-law model = 3.3x10^^ cm“^). We derive a good fit 
with xid = (1316 dof) and Nh(z) = (4.3 + 0.3) x lO^^ cm'^, 

much larger than the reported value. 


A.3. QSO B1345-h584 atz=0.646. 

QSO B1345-r584 is a Type I AGN. Corral et al. (2011) derived 
an upper limit on the intrinsic column density of Ah(z) < 1 x 10 ^^ 
cm“^ for this source. It was observed by XMM-Newton in Jun 
2001 for 40 ks, even if it was not on-axis. Because of a strong 
proton flare activity at the end of the observation only the first 
25 ks and 29 ks were retained for the pn and MOS2 instruments, 
respectively (the source falls on a CCD gap in MOSl). Data were 
extracted as above and binned to 25 counts per spectral bin. The 
fit with a power law (assuming N^^ = 1.1 x 10 ^^ cm“^) provides 
a^red - 1-19 (308 dof) with a clear under prediction of the high 
energy spectral part. In this case the intrinsic column density is 
negligible (A^h(^) <2x 10^^ cm“^). Fitting the data with a broken 
power law ,we obtain= 0.98 (306 dof). An F-test indicates 
that the addition of the spectral break is significant at > 8 cr. The 
intrinsic column density in this case is Ah(z) < 7 x 10^^ cm“^, 
in line with similar upper limits at this redshift. 


A.4. QSO B1157-1942 at z=0.450. 

QSO B1157-1942 is the lowest redshift quasar (z = 0.45) in 
the AGN-E13 sample. They quote a very low intrinsic column 
density of Nu(z) < 1 x 10^^ cm“^. XMM-Newton observed the 
quasar on Jun 2008. The exposure time is 25 ks and 24 ks for 
the MOS and pn detectors, respectively. Data were extracted as 
above and binned to 50 counts per spectral bin. Fitting the data 
with an absorbed power -law model = 3.2 x 10^^ cm“^) , 
we derive A^h(^) < 6 x 10 ^^ cm“^ (in line with previous results) 
- 1.16 (791 dof). The residuals show clear signs of an 
underlying, redshifted iron line. Adding a free Gaussian compo¬ 
nent at the quasar redshift, the fit improves to x^q^ = 1-02 (788 
dof, F-test probability > 8 cr). The derived limit on the intrin¬ 
sic column density A^h(^) < 1 x 10 ^^ cm“^ is in line with our 
expectations. 
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